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We consider a degenerate or a nearly degenerate dark matter sector where a sizable magnetic 
moment of an almost Dirac type neutral dark matter candidate N is anticipated. Then, due to soft 
photon exchange, the cross-section in direct detection of N can be enhanced at low Q 2 region. We 
discuss the implication of this type of models in view of the recent CDMS II report. 
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I. INTRODUCTION 

Recent cosmological and astrophysical observations 
strongly suggests that a significant portion of energy 
density of the universe exists in a form of dark matter 
(DM). DM has been speculated after the observation of 
the galaxy matter velocity distribution. The angular fre- 
quency distribution of the cosmic microwave background 
radiation, recent surveys of the expansion rate of remote 
galaxies and the simulation of the structure formation 
support around 25% of energy density in a form of non- 
relativistic matter while only one sixths of such matter 
has been identified as visible one. 

Among many plausible theories of DM [l[ or equiva- 
lent modification of gravity weakly interacting mas- 
sive particles (WIMP) are of particular interest [3]. The 
scenario with WIMP assumes that DM particles were 
produced in the early Universe by thermal processes and 
the DM relic density has been frozen out by decoupling 
since that epoch. The DM energy density from the cur- 
rent observation is well matched with such a thermal pro- 
duction scenario with an electroweak scale DM mass and 
appropriate electroweak scale interactions. It is also very 
interesting if DM has some interesting connection with 
the origin of the electroweak symmetry breaking as in 
the minimal supcrsymmetric standard model (MSSM) 
0. Above all, this WIMP scenario motivates searches 
for DM particles by detecting and/or producing them 
directly in the laboratory. 

Direct detection of DM has been carried out for several 
decades and the experimental sensitivity has been drasti- 
cally improved in recent years [fl Q . The abundant DM 
particles that flows through Earth may sporadically col- 
lide with ordinary matter nuclei, resulting in the recoil 
of a nucleus. The current WIMP search through the en- 
ergy deposit in the cryogenic detector device has reached 
the DM cross-section with ordinary matter nuclei at the 
order of ~ 10~ 44 cm 2 , and the upgraded CDMS II ex- 
periment 0| has reached to this level. 

In most DM models, fermionic DM is a Majorana par- 
ticle, such as neutralino of the MSSM Q. However, we 



cannot rule out a Dirac fermionic DM N at present, and 
hence it is very important to consider physical effects 
of such Dirac fermionic nature of DM. In this regard, 
we consider an appreciable magnetic moments of DM. 
In fact, the magnetic moments of neutral fermions were 
considered for a long time since the time of weak neutral 
currents [|[ up to the present age of DM [§] . 

One notable feature in the scattering through mag- 
netic moment of neutral fermion is that it has a larger 
cross section for a lower momentum photon exchange. 
The trend observed by two CDMS II candidates Q in- 
deed show this behavior: the energy deposit at 12.3 keV 
and 15.5 keV (just above the threshold of 10 keV) while 
much larger energy deposit is allowed in that experiment. 
At this time with two possible low Q 2 candidates of the 
CDMS II experiment Q, therefore, it is appropriate to 
scrutinize the Dirac DM aspect more closely. In partic- 
ular, we will pay attention to the magnetic moment / 
of A^. 1 In principle, this study includes the effects of 
the electric dipole moment also, but we will not specify 
them explicitly which would have needed an additional 
assumption about CP violation. 

There have been several works on the DM dipole mo- 
ments [T^]. In these works, various observational con- 
straints (e.g. the cosmic microwave background radi- 
ation, the cosmic 7-ray detection, the DM relic abun- 
dance) have been considered. Although the authors of 
Ref. [{| [HI commented other possibilities, they mainly 
considered DM models in which the thermal DM relic 
density is determined only by annihilation due to mag- 
netic dipole interactions. However, as pointed out by 
El , models that give rise to the DM magnetic mo- 
ment interaction with the photon exchange usually have 
other annihilation channels. Here, we consider the cases 
in which the relic density does not constrain the DM mag- 
netic dipole moment. Nonetheless, DM direct detection 
can constrain the DM magnetic dipole moment. 



f is the magnetic moment of N in units of the N Bohr magneton, 
e/2mjv ■ 



2 



As a prototype, we consider supcrsymmetic (SUSY) 
models with extra charged singlets while direct interac- 
tion between DM and nuclei can be forbidden in the 
leading order. The magnetic dipole moment is generi- 
cally induced by one-loop diagram. Such models have 
been discussed in a certain class of leptophilic scenarios 
for explaining recent cosmic ray anomalies, where DM 
particle has suppressed coupling with colored particles. 

We also emphasize that our analysis fully includes nu- 
clear anomalous magnetic moment interactions. With- 
out the photon exchange as in the neutralino case, the 
F2 form factor effects of such nucleus are usually neg- 
ligible since the range of the DM-nucleus interaction is 
much smaller than the size of nucleus. However, the in- 
teraction through photon exchange can make the contri- 
bution comparable. The nucleus magnetic moment effect 
has not been considered properly in the previous works. 
Here, we present the DM detection rates including such 
effects and compare them with the recent CDMS II data. 

This paper is organized as follows. In Sec. HH we cal- 
culate the effective magnetic moment operator from a 
simple class of supersymmetric models. In generic mod- 
els beyond the SM, one-loop induced magnetic moment 
should usually have a similar structure calculated in this 
section. In Sec. IIII| we calculate the decoupling tem- 
perature of N by solving the Boltzmann equation. We 
also use the (g— 2) e bound to constrain the hypothetical 
Yukawa coupling A. In Sec. IIV1 we calculate the direct 
detection cross-section and show the Q 2 -dependence due 
to the magnetic moment, which can be observed by the 
recoil energy distribution. The behavior due to the mag- 
netic moment is distinguishable from the predictions of 
other DM models. Next, the CDMS II experimental re- 
port is discussed in the context of DM models with a large 
magnetic moment. In this section, we also comment on 
the collider phenomenology of this class of models. Sec. 
|V] is a brief conclusion, summarizing the allowed mag- 
netic moment / of Dirac or almost-Dirac DM models. 



II. MAGNETIC MOMENT OF N 

A neutral Dirac fermion N can acquire a magnetic mo- 
ment as shown in Fig. [T] if it couples to charged par- 
ticles, fermion ip and boson 0, by a Yukawa coupling, 

rN L^+h.c. In the MSSM, the lightest supersymmet- 
ric particle (LSP) is not a Dirac fermion and thus the 
LSP DM cannot acquire a magnetic moment. However, 
in the extensions of MSSM with extra singlet chiral su- 
pcrficld, which have been proposed for many reasons : 
to address /x-problem or to raise the mass of the lightest 
Higgs boson, the LSP x °f t he MSSM sector and the ad- 
ditional singlet N may form a single Dirac fermion state 

For a specific calculation, we consider the following 
superpotential which has been discussed in (Tlj . 

W = XNe c E + m. N NN + m. E EE c + pN 3 . (1) 



In this model, there is no tree level interaction between 
N and nucleus because it couples only to leptons and its 
scalar partners. For the magnetic moment, the existence 
of the coupling of the type XNe c E is the essential one, 
which can arise in many other models extending the SM. 
For DM magnetic moment, it is required that N must be 
a DM component, presumably by an exact Zi symme- 
try. Actually Eq.(l) is the simplest model in which the 
magnetic moment can play an essential role in the DM 
detection. More complicated cases can be obtained by 
its proper extensions. 
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FIG. 1: One-loop diagrams for the magnetic moment of N. 

However, in the extension of MSSM, the interaction 
type XtpRN^+h.c. is present through NH u Hd which 
however has to be avoided by the following reasons. 
Firstly, after the Higgs fields H u and Hd develop VEVs, 
they give a VEV to N, which is harmful because of the 
mass mixing with the electron e. Second, it splits masses 
of N and N, and the anomalous magnetic moment inter- 
action is a transition type between two mass eigenstates, 
and hence the direct detection rate as a function of re- 
coil energy must be considered more carefully if the mass 
difference is of order 1-10 keV. 

For a Dirac-type neutral fermion, its electromagnetic 
interaction is dominated by the magnetic moment, as 
pointed out for the case of neutrinos Q. Let us define 
the magnetic moment of N as 

-^-Nia^NF^ . (2) 

For Eq. (JTJ, / is estimated from Fig. [U [3, 

f _ \M 2 m 2 N f 1 ^ g,(x 2 -x 3 ) 

167T 2 J I m 2 N x 2 + (to 2 - m 2 N )x + m^(l - x) 

qj,{x 2 - x 3 ) -1 

m 2 N x 2 + (m 2 , - m 2 N )x + m|(l - x) / 

(3) 



where is the mass of the fermion (e c or E) and is 
the mass of the boson (e c or E) in the one-loop diagram. 
This calculation is an illustration for a sizable magnetic 
moment of a hypothetical DM particle TV. As shown here, 
a large magnetic moment is not unreasonable since TV is 
considered to be heavy. 

If mass eigenstates splits the mass by a tiny amount, 
then the magnetic moment is of transition type with the 
initial TV and the final TV of Eq. ([2]) considered different. 
Then, the lifetime of the heavier component is 

r - 6 - bxW { f ) lloOGcvJ {Am N ) S 

(4) 

where Amjy is the mass difference between two mass 
eigenstates of this almost-Dirac fermion. If the lifetime 
falls in the 10~ 10 s region with parameters chosen appro- 
priately, then the decay products of TV deposit energy in 
the cryogenic detector. This happens for / < O(10 -4 ) in 
which case the real signal from the cryogenic data must 
be revamped. 

III. RELIC DENSITY 
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where cr ann is the TVTV annihilation cross section. 
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FIG. 2: The allowed region in the A vs. Mg plane from the 
relic density and (g — 2) M bounds for mjv = 100 GeV. The 
boundaries for mjv = 50 GeV and 150 GeV are shown as the 
lavender line and the green dash line, respectively. 



Let us now proceed to discuss the effects of a large 
magnetic moment of an extra singlet DM TV in cosmology 
and in particle phenomenology. The DM relic density is 
given in terms of the velocity averaged annihilation cross 
section (er arm v). The dominant annihilation channel of 
the Dirac DM TV allowed by Eq. (TTJ is TV + TV er +e+, 
which is mediated by exchange of the scalar component 
of E. It is straightforward to evaluate a arin v, which is 
approximately given by 



+ bv 2 +0(v 4 ). 



(5) 



a and b are 



3|A| 4 



87TTO^(1 



b = 



|A| 



B) 2 

(5B 2 - 16B - 7) 



(6) 



48Trm 2 N 



(B + l) 



m 2 g/m 2 N . Although there is also the annihi- 
lation by the magnetic dipole moment, Eq. @, we will 
neglect it because its order of magnitude is estimated as 



dipolc 
w aim u 
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N 



16tt 2 Max(mJ,mJ) 



In Fig. [2] for mjv = 100 GeV, we show the allowed 
region from the constraint on the cosmic relic density and 
also from the (<? — 2) M bound. The excluded regions are 
the lower right wedge and the upper left wedge. Toward 
the TV DM scenario, we also excluded the kinematically 
forbidden region m^r > ms- For the muon magnetic 
moment, we use the formula given in Eq. (TTJ with /i c NE 
and Eq. ([3]). The electron magnetic moment bound is 
buried in the muon magnetic moment bound. For the 
cases of w» ~ 50 GeV and 150 GeV, the boundaries 
are shown as the lavender line and the green dash line, 
respectively. 

Since / is of order |A| 2 /16tt 2 - 0.6 x 10~ 2 |A| 2 , the 
condition |/| < 10 -4 , as commented below Eq. is 
satisfied only in a small skybluc dashed-arc region in the 
lower-left allowed region of Fig. [2] Except in this small 
area, we need not consider the possibility of a heavier 
component decay in the cryogenic detector. 



IV. DIRECT DETECTION RATE 

At low Q 2 , the photon wavelength exceeds the nuclear 
size and hence the whole nucleus, with charge and mag- 
netic moment, acts as a target. Thus, the elastic scatter- 
ing cross section for the case of a spin i nucleus target 
becomes, viz. Fig. [31 
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FIG. 3: The elastic scattering of N with a nucleus through 
the magnetic moment of TV. 
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where M is the mass of the nucleus with atomic weight 
A, E Tec is the nuclear recoil energy, Q cm — Ze is the 
nuclear charge, 



A T (s,m 2 N ,M 2 



[m 2 N 



M 2 - s) 2 T Am 2 N M 2 



(8) 



and F 2 = \Ff = \ (^^f - Z) for spin § case, 

where F|* is the conventional notation. Our F 2 is a fac- 
tor i of the conventional definition, i 7 !*. The first two 
terms ( l Z terms') in the right hand side of Eq. ([7|) result 
from the electromagnetic interaction between the tensor 
operator of N (Na^ v N) in Eq. ([2]) and the vector current 
of a spin i target nucleus (ipj^ip). On the other hand, 
the last three terms ( l F 2 terms') come from the interac- 
tion between the tensor operators of N and the nucleus. 
As seen in Eq. ©, the tensor operator of a given nucleus 
defines the spin of the nucleus (tpa^ip). 

For the case of a higher spin target nucleus, '7^' and 
'cr^y' in the vector and tensor operators of the nucleus 
should be replaced by larger dimensional representations 
of the Dirac gamma matrix and the Lorentz generator, 
respectively (see e.g. Ref. [Hj]). Using Appendix, we 

i( 



can deduce that F 2 = iFf 



fj m (Z,A) / s N +l 

fi N trip 



- z 



2 1 2 2 I /ijv nip V 3s 

from Eq. (|7|) with a higher spin generalization. In the 
non-relativistic limit (£" r oc *C {M, itln} and v <C I), the 
differential cross section is given by 
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(9) 



The direct-detection rate (per unit detector mass) in a 
detector with nucleus is given by 



dR 
dEr- cc 



Pn 



d 3 vf(v) 



da 

dE la 



(10) 



Here, we assume that the WIMP of mass itin accounts 
for the local DM density and have a local velocity 
distribution f(v) with the normalization J d vf(v) = 1. 
Using a simple Maxwell-Boltzmann velocity distribution 
and pn ~ 0.3GeV/cm 3 , in Fig. 2] we plot the expected 
direct-detection rates of an almost-Dirac DM N (solid 
line) and the LSP x (dash line) in the MSSM for the spin- 
independent(SI) interactions which always dominates for 
nuclei with A > 30 in surveys of the SUSY parameter 
spaces [HI . The SI cross section is given by 

da 2 Ad , , , . , . . . 

-r^ = — (Zf P + (A-Z)f n ) 2 F 2 (E Tec ), 11 
dE rec 7TV 

where f p ~ /„ ~ 10~ 8 GeV -2 are the SI cou- 
plings of WIMPs to protons and neutrons, respec- 
tively, the SI form factor F(E rcc ) = Se^^ 32 /2 {sm(Kr) - 
kx cos(nr)) / (kt) 3 , with s = 1 fm, r = \/R 2 — 5s 2 , R = 
1.2 A 1 / 3 fm, and n = y/2m n E Icc . In Fig. IU we present 
the differential event detection rate for tun = 70 GeV, 
Eq. ([TO]), as a solid line above the CDMS II threshold 
of about 10 keV. The neutralino DM case with m x = 70 
GeV is shown as a dashed line. To show the recoil energy 
dependence, we choose an arbitrary normalization such 
that both of these lines match at 10 keV. The target is 
Ge with Z = 32 and A = 73. Its anomalous magnetic 
moment is F 2 = —36.32 as shown in Tablefl] The 1/E rcc 
dependence in the first term of Eq. ((7|) generically leads 
the dR MDM /dE lec to diverge in low E lec limit. However, 
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0(8,16) 3 


Na(11.23)i 


Si(14,28) 2 


Ge(32,73)t 


1(53,127)2 


Xe(54,131)i 


Cs(55,133)3 


W(74,183)2 




1.66812 


2.21752 


1.1218 


-0.879467 


2.81327 


0.692 


2.58 


0.117785 


F 2 


4.83 


13.36 


4.02 


-36.32 


94.56 


6.52 


83.93 


-26.30 



TABLE I: Magnetic moments /j, of several target nuclei used in cryogenic detectors. Here, fiN is the proton Bohr magneton 
jUjv = e/2m p . The nuclear spin is denoted as superscripts. F2 is the half of the conventional definition, F2 = |-Ff 1 which is 
related to the effective anomalous magnetic moment corresponding to the equivalent particle that has the same mass, charge 
and magnetic dipole moment as the target nucleus. 




10 20 30 40 50 60 70 80 90 100 

E rec [keV] 

FIG. 4: WIMP(iV)-nucleus detection rate as a function of 
the recoil energy, normalized at jf? rcc = 10 keV. The solid 
line corresponds to the magnetic moment and the dashed line 
corresponds to the SI interaction. The MDM graph and the 
SI graph meet at two points due to our normalization to their 
equality at 10 keV. So, at very low Q , the MDM is bigger 
than SI as expected. 

such a divergent effect is diminished by small velocity of 
WIMP, v 2 ~ 10 -6 which appears in the coefficient of the 
l/-E rec -term as follows: 

A_ 2M\p\ 2 n . v 2 

= ~ 2M 3 

'2.\! /'. r , < E rcc E TCC (12) 

~2M 3 . JiL- -. 

10(-B~-*)(GeV) 

Due to the suppression of the IR divergent feature 
in non-relativistic scattering, the dR MDM /dE lec diver- 
gently larger than dR° /dE Tec only in the region below 
E iec ~ lOkeV, while the slope of the dR MDM /dE iec 
above E Iec ~ lOkeV becomes eventually more flatter 
than the SI. This is why the dR MDM /dE ICC distribu- 
tion looks larger than dR SI /dE ICC above E lec lOkeV in 
Fig. (UJ. This non-relativistic v 2 - suppression in the 
l/-E rec -term might result in further interesting possibil- 
ity. Due to the suppression, 'Z terms' can be comparable 
to each other, even with 'F 2 terms'. In this regard, we 
obtain the upper-bound on the DM magnetic moment 
from the recent CDMS II data for both of the F 2 = 
and F2 7^ cases. The "maximal gap method" [l6| is 
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FIG. 5: The allowed region of the DM magnetic moment / vs. 
the DM mass mjv- The upper colored regions are excluded 
for F-2 — and F2 7^ 0, respectively, with the 90 % confidence 
level for the CDMS II data @]. 



used to estimate the proper allowed region with a 90% 
confidence level. The most stringent bound appears as 
/ < 2.88 x 1CT 3 for m N ~ 21 GeV and F 2 ^ 0. When 
we ignore the contributions of nucleus anomalous mag- 
netic moment, F 2 = 0, then the upper-bound of / is 
1.16 x 10~ 2 for Tfi 7v ~ 100 GeV. Taking into account the 
non-zero F 2 , the upper-bound goes down to 1.04 x 10~ 2 
for m^r ~ 100 GeV. The event rate becomes larger 
with non-zero F 2 so that the allowed region is more con- 
strained, producing 1 — 10% of difference in /. However, 
we can easily expect that depending on the materials in 
direct detection experiments, the difference can be sig- 
nificantly amplified due to the enhanced magnetic dipole 
moment effect. Thus, it is worthwhile to study the DM 
multi-pole interactions with nuclei more carefully. 



V. CONCLUSION 

We considered a Dirac or an almost-Dirac DM N which 
may acquire a large magnetic moment. Using the possible 
dipole interactions, we estimated the signal event rate 
which is expected in the CDMS II experiment. Using 
the recent report of the CDMS II experiment, we present 
the upper-bound of the magnetic dipole moment. The 
most stringent bound appears for ~ 21 GeV and 
F2 7^ 0: / < 1.4 x 10 -3 with a sizable nucleus anomalous 
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dipolc moment Fi contribution. We point out that this 
sizable F 2 contribution is possible in the non-relativistic 
scattering of WIMP and nucleus, leading to a several 
factor improvement in the exclusion plot. In the future 
refined direct DM search experiments, the possibility of 
DM magnetic moment of ./V can be probed with another 
independent information on its mass. 
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For the dipolc-dipole interaction contribution, the cross- 
section in the leading order becomes 



a oc tr[pM t M] 



x 



IffDMSnucleil 



(2s N + l)(2s x + 1) 



tr 



(15) 



|SDM.9nuclci| 



■sn(sn + l)s x ( s x + !) 



where the dot product is over the SO(3) space and Tr 
is over t he spin mult iplicities. Here, if we use ^nuclei = 
2M sn{sn + 1), we can see that the term propor- 
tional to /i 2 is independent of the nucleus spin. Here, if 
we use <7nuciei = 2Mfi/esN we sec that the term involving 
fx 2 is proportional to (spj + 1)/3sjv. It confirms that our 
paramctrization of Fi in the text is appropriate. 



Appendix 

In the calculation of scattering cross-section of DM 
with target nuclei, we treated the nuclei as a spin i 
particle. However, many nuclei used in the direct de- 
tection of DM have spin different from i. Moreover, 
our paramctrization of the magnetic moment using F-x, 
in which spin i for nuclei has been assumed, should be 
matched to produce a realistic value. For that, we con- 
sider two body wave function, in non-relativistic quan- 
tum mechanics with non-local interaction, which is ob- 
tained after integrating out the photon field, 

V m (x 1 y)=N a (x)® X l (y) (13) 

where N a (x) is the nuclear wave function with Sx^ = a 

and x % {y) is the DM wave function with Sz = i. The 
Scrodingcr equation for this two body system reads 

I 2M 2m N ) 2m N \x-y\ 

ffDM.9nuclcie 2 g(N) . c?(x)vj> 

2m N 2M\x-y\ 3 
= {Ha + i/ m o-di + H^i-di)^ 

where -ff mo -di and i/di-di ar e the monopole-dipole and 
dipolc-dipole interactions, respectively. Here, we ne- 
glected unimportant normalization for each terms be- 
cause we only need to match the dependency on the nu- 
cleon spin. 

Since the initial state of the problem is unpolar- 
ized, its density matrix is proportional to the iden- 
tity in (2s n + l)(2s x + l)-dimensional space, that is, 
p = 1/(2sjv + l)(2s x + 1), where sn and s x are the spin 
of nuclei and DM, respectively. Therefore, after averag- 
ing the polarizations, the cross-section is proportional to 
tT[pM*M], where M is 

M = H in t + H nlt \ r —H mt . ( 14) 
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